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MAGNETIC  COMPRESSION  OF  CHARGE -NEUTRALIZED  ION  BEAMS 
I.  INTRODUCTION 

High-power  relativistic -electron -beam  generators  have  recently  been 
used  to  create  Intense  ion  beams  characterized  by  20-50  ns,  20-200  kA 
pulses  of  .5-1.5  MeV  protons  or  deuterons  with  100-1000  cm2  cross- 
sectional  areas. 1-5  Such  beams  may  be  of  value  to  controlled  fusion  re- 
search1 in  the  areas  of  plasma  heating,6  magnetic  confinement  (by  the 
creation  of  field-reversing  ion  layers7),  and  pellet  fusion.8  Useful- 
ness in  these  areas  would  benefit  from  compression  of  the  ion  beams  to 
current  densities  higher  than  the  .1-1  kA/cm2  indicated  by  the  above 
parameters.  Two  accelerator  configurations,  the  magnetically-insulated3 
and  self-pinched4  diodes,  permit  one  to  achieve  higher  current  densities 
by  geometric  focusing.9  This  technique  cannot  easily  be  applied  in  a 
reflex-triode  geometry3  in  which  a charge-neutralized  ion  beam  is  ex- 
tracted into  a drift  region  along  a constant  axial  magnetic  field.  In 
the  work  to  be  presented  here,  the  conditions  under  which  such  a beam 
can  be  radially  compressed  by  a gradual  increase  in  the  applied  magnetic 
field  is  investigated.  A long  magnetic- fie Id  compression  length  is  de- 
sired in  order  that  the  ions  remain  closely  tied  to  the  field  lines,  i.e. 
the  laminar,  paraxial  nature  of  the  beam  flow  is  to  be  preserved  as  much 
as  practical  during  compression.  In  this  way,  the  bulk  of  the  beam 
kinetic  energy  remains  associated  with  motion  parallel  to  the  applied 


Note:  Manuscript  submitted  November  8,  1976. 


1 


I ' 


i 


field,  the  ion  gyroradius  stays  small  in  comparison  to  the  beam  radius, 
and  maximum  amplification  in  current  density  is  achieved.  In  a sense, 
the  goals  are  the  opposite  of  those  associated  with  the  creation  of  a 
rotating  ion  ring  by  passing  a cold,  paraxial  beam  extracted  from  a re- 
flex triode  through  a sharp  magnetic  cusp.7 

Cold  ion  and  electron  fluid  equations  are  solved  for  passage  through 
a specific  magnetic-mirror  configuration.  The  compressed-beam  radius  and 
the  energy  associated  with  motion  transverse  to  the  field  in  a high 
axial-field  region  beyond  the  mirror  area  are  calculated  as  a function 
of  the  mirror  ratio,  mirror  length,  applied  field  strength  and  ion-beam 
parameters.  Conditions  which  insure  the  desired  smallness  of  these  two 
quantities  are  determined.  The  effects  of  initial  beam  rotation,  con- 
vergence, and  velocity  shear  as  well  as  self-field  effects  are  considered 
in  the  analysis.  Finally,  an  example  of  how  magnetic  compression  of  ion 
beams  may  be  used  to  heat  a thermonuclear  plasma  is  presented. 
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II.  SOLUTION  OF  THE  FLUID  EQUATIONS  OF  MOTION 

The  fluid  flow  of  a cold,  charge -neutralized  ion  beam  (charge  state 
Z)  in  an  external  magnetic  mirror  field  is  now  considered.  Cylindrlcally 
symmetric,  steady-state  equations  in  the  small  electron-gyroradius  limit 
are  studied.  A self-consistent  electric  field  determined  from  the  quasi- 
neutral assumption  (n^  = Zn^  = Zn)  is  included.  Electron  temperature 
and  inertia  are  ignored.  However,  electron-temperature  effects  are  con- 
sidered in  the  appendix  where  their  neglect  in  the  present  analysis  is 
discussed.  Consideration  is  limited  to  magnetic  fields  which  change 
sufficiently  slowly  to  maintain  most  of  the  ion  kinetic  energy  in  motion 
along  the  field  lines.  This  limitation  allows  one  to  neglect  the  finite- 
gyroradius  contribution  of  the  ion-stress  tensor10  and  provides  simpli- 
fying assumptions  for  solution  of  the  equations  of  motion.  Field  and 
ion-beam  conditions  which  Insure  the  validity  of  these  approximations 
are  derived  a posteriori.  For  the  present,  the  self-magnetic  fields 
associated  with  the  beam  are  also  ignored  so  that 

V X B = V • B = 0 

determines  the  applied  magneto -static  field.  The  ion  density,  electric 
field,  ion  and  electron  fluid  velocities  (V  and  v)  are  then  determined 


from 


V.  (nV)  = V.  (Znv)  = 0 (1) 

n^V-VV  = Ze  (E  + V X B)  (2) 


3 





0 = Ze  (E  + v X B) 


(5) 


7 X E = 0 ; E = - Ttp. 


(4) 


Since  Eq.  (4)  requires  E0  = 0,  the  0-component  of  Eq.  (3), 


v B = v B 
r z z r 


(5) 


states  that  the  electron  fluid  follows  the  line  of  force.  Taking  the 
scaler  product  of  Eq.  (3)  with  B results  in 


B |$+B  |B.O 

r or  z dz 


(6) 


so  that  c p is  constant  along  a line  of  force.  Magnetic  fields  of  inter- 
est take  the  form 


-4  ■*  A 

B = 7 X (e 


so  that 


cp  = f(rA0) 


(7) 


satisfies  Eq.  (6).  The  r-component  of  Eq.  (3)  yields 


ve  = ’ VV 


(0) 


Manipulation  of  the  electron  continuity  equation,  Eq.  (3),  and 
7 X B = 0 leads  to 


-v  ^ nv 

®r  57  +Bz  h F*  = 0 
z 


so  that 
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(9) 


nv 

g(tA0). 
z 

Solutions  of  the  ion  equations  are  more  complicated.  Some  aspects 
of  analysis  are  simplified  by  noting  that  any  two  components  of  Eq.  (2) 
can  be  replaced  by  statements  of  conservation  of  canonical  angular 
momentum  and  energy  of  a fluid  element. 

(Vr  h * \ k)  ['(''o  + I7  A b]  =°  <“» 

(vr  h + \h)  (v2  + **)  - 0 

Specific  forms  for  AQ  and  cp  which  satisfy  the  requirements  of  the  pre- 
sent problem  are  chosen.  For  simplicity,  B depending  only  on  z is 

z 

chosen  so  that 


A9  = £rBz(z)  (12) 

and 

Br  = ' * rBz'  (13) 

where  the  prime  denotes  differentiation  with  respect  to  z.  The  simplest 
form  for  cp  is  then 

cp(r,z)  as  SrAg  = ^Sr^  (14.) 

where  S is  a constant. 

The  beam  is  injected  at  z = 0 into  a constant  solenoidal  field 

^ A 

B = BQez.  Solutions  are  sought  in  the  fie Id -compress ion  region  0 < z £ L 

and  in  the  high-field  solenoidal  region  z > L in  which  B = B.e  = AB  e . 

0 f z o z 
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The  plane  z = 0 may  be  Identified  with  that  of  the  transmission  cath- 
ode of  a reflex  trlode.  A radial  electric  field  In  the  z = 0 plane 
(permitted  by  Eq.  (1*0)  may  exist  when  the  cathode  is  virtual  because  of 
low-energy-electron  space-charge  accumulation  but  is  likely  to  be  neg- 
ligible when  a thin -foil  transmission  cathode  is  employed.  In  the  latter 
case,  cp  = 0 solutions  to  the  equations  of  motion  are  appropriate  since 
any  conducting  cathode  plasma  would  limit  the  radial  electric  field  to 
very  small  values. 

Assuming  a rigid-rotor  behavior  for  the  ion  flow 


VQ  = U)(z)r 


leads  to  a solution  of  Eq.  (10)  of  the  form 


V = 
r 


r 

2 w 


V 

z 


where 


w(z)  = ou  + ^ n 


and 

fl(z)  = — B (z). 
m,  z 


(15) 


(16) 


If  p(z)  represents  the  radius  of  a fluid  element  as  it  flows  along  z, 
then  integration  of  Eq.  (16)  using 


results  la 


p2!*  = constant 


Substitution  from  Eqs.  (13)  and  (14)  Into  the  z -component  of  Eq. 


(2)  leads  to 


(i?;fc  + fc)'1  = n'r2  (“-s)  • 


Employing  Eq.  (l6)  and  Integrating,  one  obtains 


V®  = V2  [1  -v(z)r®] 
z o 


where 


d_  / v\  _ _ fl^fuo-s) 
dz  \wy  ~ V^w 


and  the  constant  Vq  is  the  lon-flow  velocity  at  r = 0.  The  initial 
value  v(0)  can  be  determined  from  the  radial  electric  field  at  z = 0 
when  the  ion  motion  between  anode  and  transmission  cathode  is  purely 
axial.  In  that  case,  conservation  of  energy  requires 


i m V^rjO)  + Zecp(r,o)  = Zecp^ 

1 Z <3 


where  cp  is  the  anode  potential.  Substituting  from  Eqs.  (l4)  and  (19) 

dL 

evaluated  at  z = 0 and  collecting  terms  in  powers  of  r results  in 


V2  = 2Zecp  /m  ; v(0)  = Cl  S/V2 

o Ta  i o o 


where  0 = 0(0).  Equation  (22)  may  be  applicable  when  the  triode- 
c 

current  self  magnetic  field  is  small  compared  to  the  applied  axial 
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field.  When  motion  transverse  to  the  applied  field  is  not  ignorable, 
determination  of  V at  z = 0 requires  a complicated  analysis  of  ion  and 
electron  motion  coupled  with  a 2-dimensional  solution  of  Poisson's 
equation. 

The  ion  continuity  equation  may  be  solved  by  substitution  for 
from  Eq.  (16).  The  resulting  equation 


|-(r2nV  ) = £ (r^nV  ) 
w ozv  z 2 or  ' z 


has  as  a solution 


i(r,z)V  = K’  fCi^w) 


with  f determined  from  the  radial  distribution  of  the  incident  beam 


) = i^nCr.OV „(r,0). 
o z 


For  n(r,0)  = n = constant  when  r £ R and  equal  to  zero  otherwise,  and 
o o 

V (r.O)  = V 
z o 


n(r,z)  = 


(z)  1-v(z)j 


for  r^R(z),  while  n(r,z)  = 0 for  r > R(z).  Equation  (17)  is  the  form 


R2(z)w(z)  = R2w 
o o 


has  been  used  to  eliminate  w(z)  from  Eq.  (2?).  The  quantity  R(z)  is 

the  radius  of  the  beam  edge  with  R = R(0),  w =10  + h n , and  ol=  uj(0), 

o o o o'  " 

Determination  of  uj(z)  or  w(z)  is  all  that  is  required  to  complete 


the  solution  of  the  ion  equations  of  motion.  Manipulation  of  the 
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r-component  of  Eq.  (2)  and  substituting  from  Eq.  (18)  leads  to 

(uj2  + flu  - sft) 


n '(tu-s)r2l 


ft) 


(26) 


V 

which,  upon  evaluation  of  — from  Eq.  (16)  and  V from  Eq.  (19)  admits 

v z 

z 

of  solutions  of  the  form  <u  = uj(z)  only  when  the  cubic  terms  in  r can  be 
neglected.  The  two  terms  in  question  can  be  ignored  when 

V 

| vr2  | « 1 ; fl'r  | ^ | « 0 (27) 

z 


for  all  z s 0 and  r £ R(z).  It  will  be  shown  later  that  the  conditions 

of  Eq.  (27)  are  satisfied  for  sufficiently-slowly-varying  magnetic 

fields  provided  they  are  satisfied  at  z = 0.  The  condition  | v(0)R2  | 

« 1 is  satisfied  by  the  realistic  limitation  of  having  small  radial 

electric  fields  at  z = 0 compared  to  the  axial  accelerating  field  in  the 

triode  and  by  requiring  that  ion  motion  in  the  triode  not  be  dominated 

by  self-fields.  The  condition  I V /V  | . « 1 states  that  the  inci- 

dent  beam  is  only  weakly  convergent  or  divergent  and  is  satisfied  by 

the  same  field  restrictions  as  the  ^(OS2  condition. 

o 

Elimination  of  the  r3  terms  in  Eq.  (26)  leads  to 

(28) 

dw 

The  initial  value  of  -r-  is  determined  by  the  radial  convergence  or  di- 


vergence of  the  beam  according  to 


(29) 


since  V (r,0)  °*  V is  correct  to  lowest  order  in  Vr2.  Since  V /V  is 
z o r z 

nearly  proportional  to  r,  the  incident  beam  either  converges  to  or  di- 
verges from  a point  focus  on  the  axis  of  symmetry.  Denoting  the  axial 
location  of  this  point  by  z^,  w*(0)  is  given  by 

w'(0)  = 2w(0)/zf  . (30) 

Equations  (20)  and  (28)  may  be  solved  numerically  for  a given  set 
of  initial  conditions  and  a specified  Bz(z).  Figure  1 illustrates  the 
evolution  of  the  beam-envelope  radius  R(z)  (obtained  from  w(z)  using 
Eq.  (25)),  cu(z),  and  v(z)  for  the  case  of  parallel  (ou^  = 0;  z ^ = ®), 
shear-free  (v(0)  =0)  injection  of  ions  in  the  absence  of  electric 
fields.  The  magnetic  field  variation, also  shown  in  the  figure,  is  of 
the  form 

B 

B = [A+l  - (A-l)  cos  (TTZ/L)]  (31) 

Z d 

for  0 £ z s£  L,  while  B = AB  for  z > L.  A mirror-ratio  of  6 and  a 
* z o 

mirror  length  /V  =10  were  chosen  for  the  illustrated  calculation, 
o o 

Although  ions  are  injected  parallel  to  the  field  lines,  inertia 
causes  them  to  deviate  from  the  lines  of  force  as  the  field  is  com- 
pressed. When  the  distance  L over  which  the  field  varies  is  short,  a 
large  fraction  of  the  ion  kinetic  energy  is  transferred  to  motion 
transverse  to  the  lines  of  force  after  passage  through  the  mirror. 

Large  radial  excursions  associated  with  the  transverse  energy  reduce  the 
effective  beam  compression  since  R^,  the  maximum  of  R(z)  in  the  regime 
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z > L,  can  be  comparable  to  Rq.  When  L is  sufficiently  large,  ions  can 
closely  follow  the  gently-curving  field  lines  and  the  transverse  energy 
fraction  stays  small.  In  this  case,  R^/ Rq  approaches  A-^,  the  field-line 
radial-compression  ratio.  The  particular  case  shown  in  Fig.  1 shows  a 
5-fold  increase  in  current  density  ~(RC/RM)2  for  a 6-fold  increase  in 
magnetic-field  strength.  The  dependence  of  R^  and  the  transverse -energy 
component  on  field  and  beam  parameters  is  discussed  in  the  next  section. 
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III.  DEPENDENCE  OF  TRANSVERSE  ENERGY  AND  RADIAL  COMPRESSION  ON  FIEID 
AND  BEAM  PARAMETERS 

The  total  power  convected  by  the  beam  in  the  form  of  ion  kinetic 

energy  across  any  z = constant  plane  Is 

-RU) 

P =s  J ^m^y2*nVz*2Ttrdr 

- * "ivo  /-  [i--2  ♦ ^ * i (£)***]  • 

= 4 Vo  * + l?(jp  + i (*-)  - v]  i (32) 

where  was  evaluated  to  lowest  order  in  vr2  by  replacing  V with  V 


= / n V • 2TTrdr  ; I 
Jo  * 


= w £ n vz  * 2TTr3dr 


Substituting  from  Eq.  (23)  into  Eq.  (33) » setting  u = i^w  and  using 

Eq.  (25)  to  convert  the  upper  limit  of  integration  results  in 

/R!wn  R?w. 

f(u)  ~ ; I = TT  J f(u)  du 

o o 

where  f(u)  is  defined  in  Eq.  (23).  Therefore,  both  $ and  I are  con- 
stants which  depend  only  on  the  radial  distribution  of  the  beam  at  z = 0. 
For  the  uni form-current -density  case  described  by  Eq.  (24),  f(u)  is  pro- 
portional to  u and 

1 = * 0*0 
The  quantity  multiplying  I in  Eq.  (32)  can  also  be  shown  to  be 
constant  when  electric-potential  changes  are  negligible  compared  to  the 
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will  be  a constant  to  order  VR2.  Use  of  Eq.  (25)  then  results  in 


constant. 


(36) 


The  power  P is  then  also  a constant.  For  uniform  incident  current  den- 
sity, it  is  given  by 


P = P { 1 + 
o 


(37) 


where  Pq  = The  power  can  be  divided  into  those  associated  with 

motions  in  the  z-direction  and  perpendicular-to-z  direction  according  to 


P||(z)  = PQ  (1-fc  R2v) 


(38) 


P.  (z)  = £ P R2 
A o 


[<**(V)‘] 


* p - Pi|(z)  . 


Since  P is  constant,  Ap  = -AP(|  or 


(39) 


Vz)  = Pxo  + * P0  [r2(z>v<z>  - R2v(0)j 


where 


P,  = £ P R2 
xo  o o 


(1*0) 


(M) 
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Thus,  a knowledge  of  the  change  in  vr2  is  sufficient  to  characterize 

the  partition  of  beam  energy  between  motion  along  and  transverse  to  the 

axis  of  symmetry  of  the  beam.  In  the  region  z > L,  the  quantities  F{| 

and  P^  are  also  constant  since  VR2  is  by  virtue  of  Eqs.  (18)  and  (25). 

This  final  value  of  P is 

x 

Pi£-pxo  + *Po  fa*.  ■ vW>1?l  • <^) 


The  radial  excursion  in  the  z > L regime  can  be  determined  in  terms 
of  P^£.  From  Eqs.  (25)  and  (59) 


1 IV  + * /Vf]  _ 2P-if 
w v*  + *U j rjr?  • 

|_o  ' ' J ooo 


(43) 


Extremum  values  of  R(z)  are  obtained  by  solving  the  quadratic  equation 
in  w which  results  from  setting  w*  ~ Vr  = 0 and  replacing  m with  w-^ 

= U)  - \ AflQ.  Use  of  Eq.  (25)  then  gives 


fef'i 


l+2iu  /Cl 
o o 


(44) 


l+6T(2e-€2)?  * 

In  Eq.  (44),  is  obtained  using  the  minus  sign  while  R^,  the  minimum 
beam  radius  in  the  z > L regime,  is  determined  with  the  plus  sign.  The 
quantity  of  € is  given  by 

2 

2P,  «V  2 4p 


€ = 


^fV  z±j  (Jo  \ _1 

AP0  WV  14^, 


/Cl 
o o 


(45) 


Therefore,  the  two  primary  quantities  of  interest  P^f  and  R^R  may  be 
determined  from  initial  conditions  and  (vR£)f  which,  in  turn,  is  calcu- 
lated by  simultaneous  numerical  integration  of  Eqs.  (20)  and  (28). 


Figure  2 shows  the  limits  to  radial  excursion  and  values  of  (vR2)^ 

for  z > L as  a function  of  mirror  length  for  various  mirror  ratios. 

Initial  conditions  are  as  described  for  Fig.  1,  i.e.  v(0)  = oo^  =V^(rtO) 

= S = 0.  As  an  example  of  the  results  shown  in  this  figure,  consider 

a 10  cm-radius  beam  of  1 MeV  ions  accelerated  in  a reflex-triode  with 

B = 5 kG.  This  beam  can  be  compressed  a factor-of-9  in  area  with  less 
o 

than  a 1#  transfer  of  energy  into  transverse  motion  provided  that  the 
A = 10  magnetic  compression  occurs  over  a length  L > 4m. 

We  next  consider  how  injection  of  a beam  with  non-zero  transverse 
kinetic  energy  alters  the  compression  results  already  presented.  Figure 
3 shows  the  variation  in  RM  and  (vr2)^  with  initial  frequency  of  rota- 
tion u)  for  LH  /V  =20  and  various  field-compression  ratios  while  main- 
o o o 

taining  V^Cr.O)  = v(0)  = S = 0.  The  quantity  P^  may  be  determined  from 
the  figure  and  Eq.  03)  in  the  form 


(vR2)f(v  /Rn  )= 
Wo 


For  positive  initial  rotation  with  respect  to  the  final  value  of 
P^  is  seen  to  be  approximately  A+l  times  the  initial  value  while  R^/R^ 
is  correspondingly  larger  than  the  value  obtained  in  the  absence  of 
initial  rotation.  Increases  in  P^  due  to  initial  negative  rotation 
are  less  and  values  pf  R^/Rq  slightly  below  the  field-compression  value 
are  calculated.  It  is  also  predicted  that  values  near  - .03  re- 

sult in  final  values  of  P^  less  than  that  occurring  in  the  absence  of 
initial  rotation.  Figure  3 therefore  indicates  that  a small,  negative, 
initial  rotation  of  the  beam  (Induced  by  configuring  the  trlode  in  a way 


that  causes  a positive  radial  electric  field  to  exist  between  anode 
and  transmission  cathode)  may  allow  higher  compression  than  parallel 
beam  injection  without  increasing  the  final  value  of  P^. 

Figure  4 shows  the  variation  in  and  (vR?)f  due  to  an  initial 
radial-velocity  component  for  tu^  = S = 0.  The  initial  value  v(o)  was 
determined  by  assuming  that  Vr(r,0)  is  produced  only  by  the  self- 
magnetic  field  associated  with  the  triode  current.  If  the  radial 
electric  field  in  the  accelerating  region  is  ignored  so  that  the  plane 
z = 0 is  an  equipotential,  then  all  ions  cross  this  plane  with  the  same 
energy.  Therefore 

V^r.O)  + V^r.O)  = constant. 

IT  Z 

Substituting  from  Eqs.  (19),  (29),  and  (50)  results  in  v(0)  = z ~z 
to  lowest  order  in  v(0)r2.  The  final  value  of  P^  for  this  case  is 
then 

Pxf  - zf2  x , 

p“  = (vR2)f  r-a-  5 Zf-1  * 0 

io  Ro 

from  Eqs.  (1*0)  and  (4l).  Examination  of  Fig.  4 shows  that  for  initially- 
divergent  beams  (positive  abscissa  in  the  figure),  the  final  value  of  P 
is  approximately  AP^q  and  is  somewhat  less  for  Initially  convergent  beams 
Interestingly,  an  initial  radial  velocity  has  little  effect  on  the  final 
radial  compression  of  the  beam  for  all  but  small  values  of  A.  Some  small 
degree  of  beam  focussing  may  be  preferred  to  parallel  injection  in  order 
to  minimize  P for  high  values  of  field  compression. 
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IV.  SELF -MAGNETIC-FIELD  AND  EQUATION-OF-MOTION  SMALLNESS  REQUIREMENTS 
Conditions  are  determined  which  insure  that  the  self-magnetic  field 
in  the  drifting  beam  and  the  neglected  terms  in  the  equation  of  motion 
do  not  strongly  effect  the  calculated  compression.  The  self  field  is 
determined  by  calculating  the  total  current  density  using  results  of 
Section  II.  Denoting  the  self-field  by  b 

V x b = 4rrZen(V-v)  (k£) 


or,  in  component  form  using  Eqs.  (5),  (8),  and  (16), 


9z 


+ 2nZen*r 


Vz  “IT  vz) 


(47) 


9b  9b 

= LTTZen-r  (o»+S) 


1 9 


7 57  rb0  = ^n(Vz-v2), 


(^9) 


For  the  density  distribution  of  Eq.  (24)  and  the  condition  v (0)  = 

z 

V (°)  - V 
z o 


nV  = n V — 
z o o w 

o 


nv  = N V . 
z o o n 

o 


Substituting  into  Eq.  (47)  or  Eq.  (49)  and  integrating  leads  to 

(“T‘%)  (50> 


b8  = 2^r  I 


where  Jq  - Zen^V^.  The  long  axial  scale  lengths  compared  to  those 


for  radial  variations  allows  one  to  approximate  b^  by  neglecting  the 
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m 


wu 

term  in  Eq.  (48)  with  the  result  that 


b ®*  ^TTZenCurfSjr2  . 
z 


Then,  from  V«b  = 0, 


br  “*  ~ Ze[n(uH-S)]  'r3  . 


It  is  required  that  the  self-fields  do  not  strongly  perturb  the 
calculated  .ion-fluid  flow.  Statements  of  this  condition  from  Eq.  (2) 


VaB  » V b 
9 z 1 ’ z Q 


VJB  » V b. 
9 r 1 1 r 9 


for  all  values  of  z and  r^R(z).  Numerical  evaluation  of  velocity  and 
field  quantities  for  a wide  variety  of  initial  conditions,  beam,  and 
field  parameters  demonstrates  that  the  condition 


B n »2ttJ  V 
o o o o 


is  sufficient  to  satisfy  these  two  inequalities.  In  the  same  way, 
numerical  evaluation  of  the  quantities  in  Eq.  (27)  results  in  the 
sufficiency  condition 


3R  2A2V  « L3n  . (54 

o o o w 

Equation  (54-)  is  an  empirical  relation  derived  through  evaluation  of 
Eq.  (27)  for  various  values  of  L and  A.  Since  this  relation  insures 
| vr2  | « 1,  its  satisfaction  also  implies  small  final  values  of 
and  values  of  R^R^  close  to  A-^  provided  that  P is  small. 
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V.  DISCUSSION  OF  RESULTS 


1 


IC  has  been  demonstrated  theoretically  that  a neutralized  Ion  beam 
extracted  from  a reflex  trlode  can  be  magnetically  focussed  to  an  order- 
of -magnitude  higher  current  density  by  passage  through  a magnetic  mirror. 
If  the  incident  beam  has  only  a small  fraction  of  its  kinetic  energy  in 
motion  transverse  to  the  magnetic  field,  the  compressed  beam  will  also 
have  a small  fraction  in  transverse  motion  provided  that  the  scale-length 
L over  which  the  field  varies  is  sufficiently  large.  A general  condition 
which  L must  satisfy  is  given  by  Eq.  (54).  Final  transverse -energy 
components  for  various  initial  conditions  and  values  of  L and  A are 
illustrated  in  Figs.  2-4.  Small  final  values  of  also  insure  small 
radial  excursions  of  the  compressed  beam  about  the  radius,  Rq/A^  with 
specific  dependences  determined  by  Eq.  (44). 

A second  general  condition  which  the  field  and  beam  parameters  must 
satisfy  is  given  by  Eq.  (55).  Satisfying  this  inequality  insures  that 
the  self-magnetic  field  of  the  beam  will  be  insufficient  to  strongly 
alter  the  predicted  compression  values.  Existing  reflex-triode  experi- 
ments2 > 3 already  utilize  applied  magnetic  fields  which  satisfy  a similar 
condition  in  order  to  insure  that  the  extracted  beam  have  most  of  its 
kinetic  energy  in  axial  motion. 

An  example  of  how  the  technique  presented  in  this  work  might  be 
applied  to  a proposed  fusion-reactor  concept  is  now  presented.  Dawson, 
et.  al.11  have  suggested  that  a solenoidal  magnetic  field  of  about 
500  kG  confining  a kilometer-long,  1 cm-radlus  plasma  may  be  viable  as 
a fusion  reactor.  Ott  and  Sudan  6 argue  that  an  ion  beam  propagated 


1 


, 
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along  the  plasma  column  could  heat  the  plasma  to  thermonuclear  tempera- 
tures. They  determine  that  a 20  kA/ cm2 , 10  MeV  proton  beam  would  be 
required.  Based  on  present  experiments3  , 2 kA/cm2  could  easily  be 

extracted  from  a reflex  triode  operating  at  10  MV.  The  current  den- 
sity required  to  heat  the  plasma  and  confining  field  strength  may  then 
be  achieved  using  magnetic  compression  with  Bq  = 50  kG,  Rq  = 5.16  cm, 

Vq  = ^.36  X 10®  cm/s  (for  10  MeV  protons),  and  A = 10.  These  parameters 
marginally  satisfy  Eq.  (53)  since 


b n 

— ° Q _ k g 
2rrj  v 
o o 

Equation  (5^)  is  easily  satisfied  for  L = 3 m and  gives  Lfl  /V  “'30  for 

o o 

which  (vR2)^  » 10 "3  from  Fig.  2. 

Ott  and  Sudan  suggest  that  helical  ion  orbits  induced  by  traversing 
the  converging  field  region  may  provide  the  reduction  in  axial  energy- 
deposition  length  required  to  efficiently  heat  the  plasma.  Although  this 
motion  could  be  induced  by  decreasing  L,  radial  excursions  expected  for 
the  high  values  of  might  be  too  large  for  useful  heating  of  a 1 cm- 
radius  plasma  column.  However,  their  suggestion  of  employing  higher- 
atomic-number  ions  for  range  reduction  would  allow  one  to  magnetically 
compress  the  beam  in  a way  which  maintains  most  of  its  energy  in  axial 
motion  and  limits  radial  excursions  to  values  small  compared  to  1 cm. 
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VII.  APPENDIX 

Here,  the  inclusion  of  electron  pressure  terms  in  the  equations  of 
motion  are  considered.  The  electron  force-balance  equation  then  takes 
the  form10 

«■ • 

Zen(E  + v X B)  + 7»"3  = 0 (Al) 


where 


Pj.  1 + (pn-Pi>  %% 


(s)  = ™ (^-)- 0 


Here,  ^ = B/B  is  a unit  vector  in  the  magnetic-field  direction.  Since 
(V-H)Q  = 0,  Eq.  (5),  which  states  that  the  electron  fluid  follows  the 
lines  of  force,  still  applies.  Taking  the  scaler  product  of  Eq.  (Al) 
with  B yields  after  some  algebra  and  generous  application  of  Maxwell's 


equations 


(Brt?  + B*fc)  v-(\k+\h) 

+ Vi(Bri?+\  s)  B = 0- 


Employing  Eq.  (A3)  allows  this  to  be  reduced  to 

(Br  k + Bz  k)(r  - z°>)  * I (\k  + \ k)  r- 0 

-»  _ p phBs 

since  B*7t  = 0 follows  from  v«7t  = 0 and  Eq.  (5).  Since  — and  — g- 

nB  n 
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are  constant  along  a line  of  force,  and  n is  roughly  proportional  to 
B from  Eq.  (25),  then  p^  ~ B2  (®*  B2  for  the  case  at  hand)  while  pu~  B 
as  the  beam  is  compressed.  If  electrons  in  the  injected  beam  are 
characterized  by  p « p^  at  z = 0,  then  p^  will  greatly  exceed  p|(  in 
the  compressed  beam.  Also,  since  some  tens  of  electron  volts  of  temp- 
erature might  be  reasonable  for  the  injected  beam,  ^ m^V^2  >:>  at 

z = 0.  Thus,  if  important  at  all,  electron  pressure  effects  will  only 
be  significant  in  the  compressed  beam  where  px  >>  p|(. 

With  the  considerations  of  the  last  paragraph,  it  is  reasonable  to 
neglect  p in  comparison  with  p^  in  Eq.  (A5)  so  that 


or 


(A  6) 


Zecp 

n 


g — B 

2 z 


(A7) 


in  analogy  with  Eqs.  (6),  (7)  and  (14). 
so  that 


But,  from  Eq.  (A3),  — ~ B 


Zecp 


+ kT 

xo 


B 

B 


(A3) 


where  a constant  temperature  profile  at  z = 0 has  been  chosen. 

It  is  seen  that  the  maximum  change  in  Zecp  associated  with  initial 

electron  temperature  is  of  order  A*kT  which  is  much  too  small  to 

xo 

effect  ion  motion.  The  electron  motion  is  only  effected  in  azimuth 
through  the  temperature-dependence  of  the  radial  electric  field  which 
impacts  only  on  the  small  self-field  of  the  compressed  beam.  It  is 


2k 


Fig.  3 — Variation  of  final  maximum  beam  radius  and  velocity-shear 
parameter  with  initial  rotational  frequency 


-.0.2  -0.1  0 0.1  0.2 


(V0/rn0)(Vr/V2)0  — 

Fig.  4 — Variation  of  final  maximum  beam  radius  and  velocity-shear 
parameter  with  initial  beam  convergence  or  divergence.  Positive  values 
of  the  abscissa  indicate  divergence. 
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